Broken discrete symmetries in a frustrated honeycomb antiferromagnet 
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We study the magnetic phase diagram of the J\-Ji Heisenberg antiferromagnet on a honeycomb 
lattice at the strongly frustrated point J2/J1 = 1/2 using large-scale Monte Carlo simulations. 
At low temperatures we find three different field regimes, each characterized by different broken 
discrete symmetries. In low magnetic fields up to h c \/J\ ~ 2.9 the Z3 rotational lattice symmetry 
is spontaneously broken while a 1/2-magnetization plateau is stabilized around h C 2/Ji = 4. The 
collinear plateau state and the coplanar state in higher fields break the Z4 translational symmetry 
and correspond to triple-g magnetic structures. The intermediate phase h c i < h < h C 2 has an 
interesting symmetry structure, breaking simultaneously the Z3 and Z4 symmetries. At much lower 
temperatures the spatial broken discrete symmetries coexist with the quasi long-range order of the 
transverse spin components. 

PACS numbers: 



The search for a quantum spin liquid - an insulating 
magnet with a gapless ground state which breaks neither 
lattice nor spin symmetries - has been the focus of many 
studies on two-dimensional frustrated quantum antifcrro- 
magnets [IHH . Such systems are assumed to be the main 
candidates to describe a rich variety of unconventional 
phases, phase transitions and critical points with decon- 
fined fractional excitations @-|z[. Frustration plays an 
important role in classical systems as well. Within this 
context, the phenomenon of order by disorder (OBD) Q 
is the perfect example where the interplay of frustration 
and fluctuations produces the emergence of unexpected 
order. OBD implies that certain low-temperature spin 
configurations are favored by higher entropy rather than 
by lower energy. For instance, some frustrated spin mod- 
els may exhibit magnetization plateaus even at the clas- 
sical level 0, [l(| . In this case fluctuations are responsi- 
ble for stabilizing particular collinear spin configurations 
that have softer excitation spectra compared to a gen- 
eral noncollinear spin state. Another nontrivial fluctua- 
tion effect is a finite-temperature transition in 2D Heisen- 
berg antiferromagnets related to breaking (discrete) lat- 
tice symmetries in the absence of a long-range magnetic 
order [ll[ . Studying various fluctuation- induced types of 
magnetic order is important in order to establish robust- 
ness of a hypothetical spin-liquid state. 

Here, we consider the Heisenberg antiferromagnet on 
a honeycomb lattice, model realized in a number of real 
magnetic materials [ 1 2|— 1 1 51] . Recent interest in this model 
is largely motivated by the experimental realization of 
the spin-liquid state in Bi 3 Mn40i2(N0 3 ) [II]. Mulder 
et al. [ItJ have studied the J1-J2 frustrated honeycomb 
antiferromagnet, Fig. [TJa), in zero magnetic field. They 
have found that quantum fluctuations select a family of 
special states characterized by three inequivalent Q vec- 
tors for 1/6 < J2/J1 < 1/2 and by three different ones 



for 1/2 < J2/J1. While at finite temperatures thermal 
fluctuations melt the spiral order, the discrete Z3 lattice 
rotational symmetry is still broken at low temperatures. 

In this Letter, we extend the previous theoretical work 
on the J1-J2 honeycomb antiferromagnet to finite mag- 
netic fields. This is important in view of the experimen- 
tally observed field-induced transition between the spin- 
liquid state and a lon g-ra nge ordered magnetic structure 
in Bi 3 Mn40i2(N03) [l6(. Specifically, we focus on the 
highly frustrated point at J2/J1 = 1/2. Surprisingly for 
such a simple model, we find a plethora of new emergent 
broken symmetries, which may exist in magnetic fields 
alongside with a 1/2-magnetization plateau. 

Our main results are summarized in Fig. [ljc). The 
magnetic phase diagram at low temperatures is divided 
into three regions. At low fields, and very low tempera- 
tures, the spins show a quasi-long-range order (QLRO) 
in the XY plane in a canted antifcrromgnctic (single- 
q) structure which breaks the translational and orienta- 
tional symmetries of the lattice. By increasing the tem- 
perature while keeping fixed the value of the magnetic 
field, we expect a Kosterlitz-Thouless (KT) transition to 
a phase in which the internal structure melts down. This 
intermediate phase has the sp atial symmetry broken, 
similar to what was found in [171 ]. 

When the value of the magnetic field exceeds h Cl / J\ ~ 
2.9, we find a first order transition into an intermedi- 
ate phase similar to the 2-1-1 phase observed in the 
pyrochlore antiferromagnet ZnCr204 [3] which trans- 
forms at h C2 /Ji — 4 into the collinear 3-1 structure 
via a second-order transition. Although in all high field 
phases the unit cell of the system have 8 spins, we shall 
use the same terminology as in [lj|] to emphasize the 
magnetic structure for each frustrated four-spin block. 
The collinear 3-1 state corresponds to the magnetization 
plateau with M/Af sat = 1/2. It preserves the XY rota- 



2 




Canted 2-2 



2-1-1 QLRO 



3-1 QLRO 



h/Ji 



collinear phase 




2 /i ri ft-4fe+ 



FIG. 2: (color online) Magnetization curve, absolute value of 
total S z and susceptibility as a function of the magnetic field 
for J2/J1 = 1/2 for two system sizes (N = 2 * L 2 ). Inset: 
& spin configuration, uuuuuudd state (collinear 3 — 1), on 
the M/M sa t = 1/2 quasi-plateau at T/Ji = 0.01. The blue 
and green dots correspond to a spin fully polarized along the 
magnetic field and in the opposite direction respectively. 



FIG. 1: (color online) (a) Honeycomb lattice. Blue (Red) 
line represents first (second) nearest neighbors, (b) Elemen- 
tary tetrahedrons used to write the Hamiltonian in the highly 
frustrated point, (c) Schematic phase diagram. The dotted 
line represent a constant low temperature line which can be 
obtained by increasing the magnetic field. By moving along 
this line one should observe phases I, II and III with broken 
Zz, Zz x Zi and Z4 symmetries respectively. In phase (III) 
the lattice orientational symmetry is restored. If the temper- 
ature is low enough (as depicted in the figure) QLRO phases 
should be present everywhere except in the vicinity of h — 
and h C2 / J\ = 4 where the pseudo plateau with a collinear con- 
figuration is present. Finally, increasing further the magnetic 
field should drive the system into the trivial paramagnetic 
phase characteristic of the high temperature region. 

tional symmetry about the field direction but breaks in a 
special way the discrete Z 4 lattice translational symme- 
try (see the inset in Fig. [2]). For this value of the magneti- 
zation (magnetic field) increasing temperature produces 
a phase transition into the paramagnetic phase via a con- 
tinuous transition in the universality class of the Z4 clock 
model. It is interesting to notice that the 2-1-1 state has 
a supplementary Z3 broken symmetry with respect to the 
3-1 state. One would then expect again a 3-states Potts 
transition by keeping the temperature fixed and increas- 
ing the magnetic field in order to pass from the 2-1-1 to 
the 3-1 state. 

By increasing the applied magnetic field beyond the 
1/2 plateau, we enter into the high-field phase where a 
3-1 state known previously for the pyrochlore antifcrro- 
magnet [l9| becomes stable. The classical configurations 
for the spins are now canted again, we have broken trans- 
lational symmetry due to the long range order in the z 
component of the spins and we recover quasi-long-range 



order for the XY components. Note that, in contrast to 
the low field phase, the lattice orientational symmetry is 
unbroken. 

To see all this in detail, let us now introduce the model. 
The spin Hamiltonian is given by 

H = 3 X £ 8,-Sj + Ja £ S.-S, h^s; (1) 

<i,j>i <i,j>2 i 

where < i,j >i and < i,j >2 denote nearest and next 
nearest neighbors respectively, J\, J2 > and h is the 
external magnetic field. 

Let us briefly discuss the magnetically ordered phases 
of the model {1} . For small diagonal exchange J2/J1 < 
1/6, classical spins form the Neel state. For 1/6 < 
J2/J1 < 1/2 the classical ground states are degenerate 
spirals forming a closed contour around T-point, while 
for 1/2 < J2/J1, the closed contours are centered around 
K and K'-points [II]. 

We now focus on the highly frustrated point J2/J1 = 
1/2, where the Hamiltonian can be written, up to a con- 
stant term, as a sum over elementary tetrahedrons that 
we label as V and A respectively (see Fig.[]J>): 

o=v,a v 7 

where S v = Siev S «' S ^ = J2ie/\ S « and h = By 
minimizing the energy on each A and V, one obtains the 
constraint: S v = S A = h/(2Ji). The classical ground 
state is obtained when this constraint is satisfied in every 
block and presents only the typical global rotation as a 
degeneracy. The saturation field h s is determined by the 
condition = 4 and = 4 which gives h s = 8 J\. At 
this value all the spins are aligned with the z— axis. 
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Monte Carlo simulations have been performed using 
the standard Metropolis algorithm in combination with 
the microcanonical over-relaxation steps, see [HI for fur- 
ther details. Periodic boundary conditions were imple- 
mented for N = 2 x L 2 site clusters with L = 24-72. 
At every magnetic field or temperature we discarded 10 5 
Monte Carlo steps (MCS) for initial relaxation and data 
were collected during subsequent 2 • 10 5 MCS. The error 
bars were estimated from 20 independent runs initialized 
by different random numbers. 

Let us now discuss various physical quantities used to 
clarify different phases and corresponding transitions. In 
the first place, we calculate the magnetization, suscepti- 
bility and absolute value of S z defined as 



1 ^ „„ dM 

i—l 



\s z 



1 N 



(2) 



In Fig. [2] we show the magnetization curve, susceptibility 
and absolute value of S z as a function of the external 
field at temperature T = 0.01 J\. 

The susceptibility Xm shows a dip around h = h C2 = 
4 Ji which indicates the presence of a quasi-plateau 
phase. In the same region the absolute value of \S Z \, 
which measures how "collincar" is the magnetic configu- 
ration, is close to one and therefore the magnetic phase 
established is a "collinear phase" . The situation is com- 
pletely different for small fields, \S Z \ is smaller that 1/2 
and then the phase corresponds to a canted 2-2 antifcrro- 
magnetic (AF). Both regions are separated by a big jump 
in \S Z \ around h = h Cl ~ 2.9 Ji indicating a first order 
phase transition. 

The previous results suggest that the low field phase is 
continuously connected with the zero field case studied 
Fluctuations select a commensurate wave vector 



17] 



corresponding to the M-point in the Brillouin zone (BZ). 
It has residual triple degeneracy. At zero magnetic field 
the selected structure is described by a single wavevector. 
To detect this single-g-paramagnetic phase transition we 
introduce a local complex order parameter A Q/3 (r) [ItJ 
and its averages as 



A"* 3 (r) = l -S<X (r) [Sg (r)+ W Sg(r+b)W S P E 
— J2 A ra (r) + A»(r) + A zz (r) 



(r-a+b) 



|_L^A-(r)+A^(r) 



(3) 
(4) 
(5) 



where the sum over r runs over one of the two sub- 
lattices (say A), uj = e 27 ™/ 3 , a = x,y,z and a, b 
arc the primitive translation vectors of the direct 
lattice (See Fig. [T|). In the previous definition we 
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FIG. 3: Binder cumulant associated with the order parame- 
ter A as a function of temperature T showing the transition 
between the paramagnetic phase and the single-g phase at 
h/Ji — (Top) and h/J\ = l (bottom), for system sizes up 
to L = 72. 



have normalized to 1 the case of a perfect "two- up, two- 
down" collinear configuration (only achievable at T = 0). 

Following the standard procedure, the second-order 
transition between a paramagnetic phase (large-T) and a 
single-q phase may be located by the crossing point of the 
corresponding Binder cumulant Ua measured for differ- 
ent clusters. We have used instantaneous values of ([3]) to 
measure the susceptibility xa and the Binder cumulant 
Ua associated with this order parameter defined as 



XA 



Nc 
T 



<(A) g 



((A) 4 ) 
((A) 2 } 



(6) 



where N c is the number of unit cells. We illustrate this 
method in the top panel of Fig. [3] for the transition be- 
tween the paramagnetic state and the single-g state. The 
alternative approach is to study the susceptibility since 
the critical exponent rj is known precisely, r\ = 4/15 (20| . 
In the critical region the susceptibility scales as 



XA = L 2 ^f(\r\L 



= 1 - T/T c 



(7) 



Hence, the normalized susceptibility xa/L 2 ~ v becomes 
size-independent at r = and curves for different L plot- 
ted as functions of T exhibit a crossing point, similar to 
the behavior observed for the Binder cumulant. 

In the case of non zero field we find three different re- 
gions, as it is schematically depicted in Fig. JTJ): 
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(i) For h < h Cl we have a similar situation as for zero 
field, namely at fixed magnetic field and coming from 
high temperature, a paramagnetic- Z3 transition occurs. 
Decreasing further the temperature one then should en- 
counter a KT transition to a canted 2-2 AF single-g 
QLRO in the XY spin components. 

(ii) For h Cl < h < h C2 we have a finite-temperature tran- 
sition related with the breaking of Z4 x Z3 spatial sym- 
metry. Again, at low temperatures one should find a 
KT transition to a QLRO phase in the XY spin plane 
corresponding to a coplanar configuration similar to that 
found in [lj| for the pyrochlorc lattice, which we dubbed 
the 2-1-1 phase. Note that close to h Cl , thermal fluctu- 
ations select a collincar configuration for the spins which 
dramatically decrease the effective spin stiffness in the 
XY plane when entering in the phase II. Then one ex- 
pects the KT transition temperature in region II to be 
lower than the one of region I, as schematically depicted 
in Fig. p. 

(iii) For h > h C2 , coming from the high temperature 
phase, we encounter a continuous transition to a phase 
with broken Z4 translational symmetry, see the inset in 
Fig. [2 This 3-1 phase as well as the 2-1-1 phase are 
characterized by mixing of three wave-vectors (triple-g 
structures). Again, at even lower temperatures, we ex- 
pect a KT transition to a QLRO state where the spins 
adopt a planar configuration in which three spins are 
pointing in the same direction and the XY component 
of the remaining spin compensates the sum of the three 
first ones. 

Thermal fluctuations have a strong effect in fields 
around hh s , where they stabilize a collinear 3-1 state 
and there is a symmetry breaking related to this selec- 
tion, as we explain now. At this point the spin pattern 
consists of 8 spins per unit cell (Fig. [2|). We rewrite the 
coordinates of the 8 spins as we show in the inset of Fig. 
2] and introduce the following Z4 order parameter 



I 4 

nzZzZ' 

Z7 n=l 



r(n—l)gz 



r(4-n)g 



(8) 



Using this order parameter one can construct the corre- 
sponding Binder cumulant Uz 4 in the usual way. The 
results are shown in Fig. [4] (at h/J\ = 4.2), measured 
for different cluster sizes at the transition between the 
paramagnetic and the Z± symmetry breaking collinear 
3 — 1 state. Results for various temperatures and mag- 
netic field scans arc summarized in the phase diagram 
presented in Fig. QJc) 

To summarize, we have studied the phase diagram of a 
strongly frustrated classical J\-Ji Heisenberg antiferro- 
magnet on a honeycomb lattice in a magnetic field. We 
have found a very rich low temperature phase diagram 
showing three non-trivial regions characterized by differ- 
ent broken lattice symmetries, as summarized in Fig. (1). 

In order of increasing magnetic field at fixed (low) tem- 
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FIG. 4: (color online) Binder Cumulant corresponding to the 
Z4 order parameter at h/Ji = 4.2 defined in the main text. 
The transition between the paramagnetic phase and the 3—1 
state is clearly observed. 



perature, one first encounters a first order transition trig- 
gered from the low field phase (I), which breaks the ori- 
cntational Z3 symmetry, into the intermediate field phase 
(II) where a supplementary Z4 symmetry, related to lat- 
tice translations, is broken. By increasing the field fur- 
ther, a continuous transition to phase (III) occurs, which 
provides an example of (lattice) symmetry restoration. 

All these discrete broken symmetry phases should co- 
exist with the corresponding QLRO phases, after a KT 
transition occurs at lower temperatures. 

We have explicitly numerically checked that the tran- 
sition from phase (III) to the paramagnetic phase is con- 
tinuous and one should then in principle expect continu- 
ously varying exponents associated with the universality 



class of the Z4 Potts model [2l|. The richness of this 
phase diagram illustrates the importance of frustrated 
and competing interactions and the onset of fluctuation 
mechanisms in the selection of the low energy configura- 
tions. 

The present study may be relevant in the study of dif- 
ferent compounds that are described by the frustrated 
hexagonal Heisenberg model, such as Bi3Mn40i2(N03) 
fll| already mentioned and the family of compounds 
BaM 2 (X0 4 ) 2 with M = Co, Ni and X = P, As, which 
consist of magnetic ions M arranged in weakly coupled 
frustrated honeycomb lattices with spin S = 1/2 for Co 



and S = 1 for Ni [22[. In the case of materials with 



spins higher than 1/2, magnetic field experiments could 
unravel some of the structures found in the present pa- 
per. Last but not least, the controlled setup of optical 
lattices for cold atoms would allow to create arbit rary 
lattice structures as well as to tune the interactions [23J. 
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